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jective and surjective onto Z? Such a function, which we call harmonic 
labeling, is constructed when the graph is the Z 2 square grid. It is shown 
that for any finite graph G containing at least one edge, there is no har- 
monic labeling of G x Z. 



Q '■ 1 Introduction 

U' 



Let G = (V, E) be an infinite graph of bounded degree. 
Definition 1.1. A junction <j> : V(G) — > Z is harmonic if 



(N 

// 4> is also injective and surjective we say that <fi is a harmonic labeling oj 

p: G . 

In this note we study the somewhat ad-hoc question, which graphs admit 
harmonic labeling? In view of some of the techniques used, we hope that a nice 
^O | structure will emerge, making this project natural and interesting. The study 

linking between the behavior of harmonic functions on spaces and geometric 
properties of the spaces, started with Liouville theorem, is well developed, see 
e.g. [YS94] . Here we want to suggest a variant in the spirit of additive combi- 
natorics over graphs, see e.g. |AABL09] , 



In the coming sections we will prove some positive and negative results re- 
garding the existence of harmonic labeling. We are still far from classifying 
Cayley graphs or vertex transitive graphs, with respect to this property or un- 
derstanding the structure of such functions. In particular we don't know if 
admitting harmonic labeling is a quasi isometric invariance for Cayley graphs? 
It will be of interest to verify on natural examples of Cayley graphs whether 
they admit harmonic labeling and come up with useful invariants. In the last 
section we list few open problems. Before moving to trees, grids and products, 
here are two quick warm up examples. 

Example 1. The easiest example oj harmonic labeling is when V(G) = Z and 
E{G) = {{i,j) '■ \i — j\ = 1}. Let (f> = Id. Each vertex is oj degree 2, and 
i = <fi(i) = \{(f>{i + 1) + 4>{i — 1)) = \{i + 1 + i — 1) = i- For a graphical 
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Figure 1.2: No labeling on the cross 



representation of the H-labeling see Figure ["PI In the graphical representation 
we write the image of <j> on the vertices of the graphs. 

"~© © ® © 



Figure 1.1: Labeling on Z 



The next example shows that the existence of harmonic labeling is not trivial. 



Example 2. Let G be the infinite cross, as in Figure HOI If (j) is a harmonic 
labeling then so are ip(x) — (f>(x) —(f)(0) and ip( x ) = — <t>{ x )- Thus we can always 
assume zero value on some appointed vertex. Assume by contradiction that a 
harmonic labeling of the cross exists, and assume zero value at the intersection. 
From the harmonic property three values a, b, c span the value of <f> on the cross. 
We can assume that two of these are positive (if not look at ip(x) = —<j)(x)). 
Without loss of generality assume a and b are positive. We will obtain the 
number a ■ b in both the left and the right arms of the cross contradicting the 
assumption <j> is injective. 

From the examples and the the study below, the notion of label spanning 
sets (Remark 2J Section [S]) emerges as a key concept. It is of interest to study 
and classify these sets. 




Figure 2.1: Four- regular tree 



2 Trees 



Theorem 2.1. Let G be a d-regular tree, then there exists a labeling on G. 



Proof. We start with the four-regular tree. Observe Figure [2TTI 

We build : V(G) — > Z recursively. Assume <fi is defined on some subset 
V C V(G). Let 

m = min{|z| : z (jL 4>(V')}. 

, if m £ (j>(V') take — m. Define <f> to be m on one of the free leaves on the 
boundary of V. The parent of this vertex has two other children, denote by 
/ = 0(parent), / = ^(parent of parent). Those children give us the required 
degree of freedom to make <j> injective. Define <j> on one of them to be a number 
a much larger than 

max|0(w)|. 

Thus all that is left is to show that the value of the second child b is different 
than all the previous values of <fi. 
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It is easy to see that if a is large enough b is small enough and <f> is injective. 
Subjectivity is obvious from construction. 

Notice that for d > 4 we will have more degrees of freedom, and the recursive 
construction will work as well. For the three-regular tree the previous construc- 
tion fails since there are not enough degrees of freedom in each level of the tree. 
Thus the construction will use two levels in each step. Observe Figure |2~121 We 
build <p recursively, very similar to the d > 3 case. Assume <f> is defined for 
some subset V C V(G). We wish to define </> on one of the leaves with minimal 
distance to the root. Let m be as in the previous construction, a number with 
minimal absolute value that is not in 4>(V). Define to be 6 in the desired leaf, 
m on one of its children and b on the second child. Let / = </>(parent), then 




Figure 2.2: Three-regular tree 



b = -|(a + / + to). If a is much larger than all the values in 4>(V), b will be 
much smaller and </> stays injective. Once again <f> is surjective by construction. 

□ 

Theorem 2.2. Every tree such that all vertices have degree greater than 2 has 
a harmonic labeling. 

Proof. Choose a vertex and denote it as the origin F(0). Let 

r(n) = {r(n)^ : i = 1, . . . , \T(n — l)[,j = 1, ... , degree of common father — 1} 

be the set of vertices of distance n from the origin. We build the harmonic 
labeling recursively. Suppose we defined 4> on T(n — 1) lets expand it to T(n). 



Let I be the degree of common father of r(n)j, and let toi, 



,mi- 2 



be the 



numbers in N with smallest absolute value such that mi, . . . , to/_2 ^ U" =0 T(i). 
Define T(n)j = nij for j = 1, ...,/ — 2, now define F(ri) 1 _1 be big enough such 
that both it and T(n)[ be different than Uf =o r(j)\r(n)'f 1 ur(n)' 1 . We continue 
the same way for Tl, i ^ 1. 

□ 



3 Finite graphs 

Theorem 3.1. Let G be a finite graph with at least one edge. The graph G x 
does not admit harmonic labeling. 




Figure 3.1: Ladder graph Z x Z 2 

As an example to this theorem we will show that there is no harmonic 
labeling for the ladder graph. 

Proposition 3.2. The graph Z x Z2 does not admit harmonic labeling. 

Proof. Assume by contradiction that there exists a harmonic labeling 



In Figure [3~T1 we denote by 



{bijiez 



m x {0})} i& 
{<j>{i x {l})} ie] 



(3.1) 



By the harmonic property of the labeling, (f> is determined by its values ao, a,\, bo, 01 • 
Without loss of generality assume a\ = maxjao, a\, bo, bi}. Since we assumed 
that </> is injective, 61 and ao are smaller than a\ by at least 1 and 2. This yields 
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3ai — b\ — ap > 3ai — {a\ — 1) — {a\ — 2) = a\ + 3. 



(3.2) 



This means that either a 2 or 02 are the maximum of {ai, &i}f =0 and greater than 
a\ + 3. By induction we get that for every k £ IN \ {0} cither ak or bk are the 
maximum of {ai,bi}^ =0 and greater than ma.x{ai,bi} i ~Q + 3. Assume that at 
least one of the next limits are true liminf n _ ) . 00 a n = —00, liminfn-yoo b n = —00. 
Without loss of generality assume the second limit is true. There is some Ni € N 
such that bjsix = minjai, bi} i= t Q . Now for all fc 6 N, bjsix+k < ^JVi+(fc-i) — 3 and 
&Ni+k > awi+(fc-i) + 3- If non of the above limits are true the the right side of 
the ladder is bounded from below and this leaves us with two options. 

1. The maximum value changes infinitely many times between a and 6. 

2. From some i, bi is the minimum of a,;, bi, aj+i, 6j+i, and thus for each j > i 
bj < bj+i — 3. 

In any case The density of positive integers in the right part of the ladder is less 
than 2/3. 

If without loss of generality limsup i _ > _ 00 a, = 00, there exists some N2 E — N 
such that ajv 2 = max{a.;, bi} i J:_ 1 and biy 2 = min{ai,bi} i J : _ 1 . As before for all 
fe E IN, bN x ~k < OiVi-(fc-i) — 3 and aN 1 ~k > liVi-(fc-i) + 3. Thus the density of 
negative integers is no more than 2/3. This contradicts our assumption that <p 
is surjective. □ 




c=(L+2I)b-a 



Figure 3.2: Finite graph G times Z 



Proof of Theorem \3. 11 By the harmonic property, determining <f> on two adja- 
cent slices of G determines the function on all G x Z. For a graph G with n 
vertices define the Laplacian matrix L = (lij) n x.n 




if i = j 
if i ~ j 

otherwise 



(3.3) 



L is a real symmetric matrix, and hence has a real orthogonal diagonalization. 
Let fj-o < fJ-i < • ■ • < Hn—i be the eigenvalues of L. It is known that Vi, /Uj > 
0, no = (the vector (1,1,..., 1) is an eigenvector with eigenvalue 0). The 
multiplicity of the eigenvalue is the number of connected components in G, 
and then for every connected component there is the corresponding 0-cigcnvector 
which is the vector all of whose entries are 0, accept for those which stand for 
vertices in that connected component - and their value is 1. 

Denote by y = {y;}™ =1 the values of <j) on some slice of G and x = {xi}™ =1 
the values on the slice above it. 

Lemma 3.3. Leta,b,c be the values of (j) on three slices of G one on top of the 
other. Then 

' c\ = f L + 2I -I \ / b 

b { I Ma 



(3.4) 



Proof. Immediate from the harmonic property, note that we add 21 because the 
degree of each vertex in G increases by 2 after multiplying by Z. See Figure [ 
for clarification. 



Next we analyze the eigenvalues of A = 



L + 2I 
I 



an eigenvalue of A with an eigenvector ( _ ) . Then 
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□ 
Assume A is 



Lx = | A + - - 

Thus x is an eigenvector of L with eigenvalue [i = A + j — 2. On the other 

a \ , ( a 



direction, if /i is an eigenvalue of L with eigenvector a then . _ and i 

V Xa J \ \ a 

are eigenvectors of A with eigenvalues A and i-, that satisfy /j, = A+ -r — 2. We 

know /x > so A, i = i ((fi + 2) ± ^/(/i + 2) 2 - 4) > 0. 

In particular, for \x ^ 1, there are two different eigenvalues A, t 7^ 1. 

For ii = 0, A = 1, the corresponding eigenvector is _ I. where TJ is the 

eigenvector of L corresponding to the eigenvalue and the relevant connected 

component. 

It is easy to verify that for any such v, we have that 

'(;)-(r)-(s) + C) 

Thus, we have found a Jordan decomposition for A. For any eigenvalue 
of A, different from 1, we know its eigenvector, and for A — 1 we have found 
the eigenvector and the extended eigenvector (there might be several, due to 
multiplicity). The above vectors, together, form a basis for the vector space, 
which is of dimension 2n, where n is the number of vertices of G. 

Let {vi} be the eigenvectors and extended eigenvectors of A with eigenvalues 
A;. Any vector of 2n coordinates can be represented as a linear combination of 
these vectors, and in particular any vector which represents the values of (j) on 
two consecutive slices. 

Lemma 3.4. The representation of vectors which represent the values of <fi on 
two consecutive slices, using the base {v^, uses only eigenvectors or extended 
eigenvectors with eigenvalue A = 1. 



Assuming Lemma 13.41 we can finish the proof of the theorem. Indeed, the 
basis elements which may be used for the representation of such vectors, must be 
constant on every set of coordinates which corresponds to a connected compo- 
nent of G (note that for any such component, there are two sets of coordinates, 
among the In coordinates, that correspond to that component - as there are 
two slices). This is due to the fact that both the eigenvectors of 1 (in A) and 
the extended (Jordan) eigenvectors are constant on such sets. Thus, any com- 
bination of such vectors will be constant on connected components (in any slice 
separately), but since there is at least one edge in G, there is at least one such 
component with at least two vertices. <j) must receive the same value for both 
vertices - in contradiction with <f> being bijective. 

We now continue to the proof of the lemma. 



Proof of Lemma \3-4\ It is enough to prove there are no eigenvectors for eigen- 



values larger then 1 in the representation of such a vector I _ j as inverting 

the order of the slices takes an eigenvalue A to j and vice versa and leaves 
us with another harmonic labeling. To see this remember equation (13.51) and 
calculate 



Ly + 2y — x 
V 

jx J A V x 



(3.7) 



We are left with showing that if u = I _ I represents the values of </> in two 

consecutive slices, then in it's linear combination with respect to the basis vl, 
no eigenvector for eigenvalue greater than 1 appears. We define a line to be a 
set of vertices in G x Z with the same G-coordinate. We have n lines. 

Definition 3.1. Let S = {Sn}™=?oo be a monotone increasing sequence of 
(different) integers. Denote by [n] the set {— n, ...,n— l,n}. If 



lim \[m}C]s\/(2m + l) 



m—too 



exists, it is called the density of S . 

Some well known (and easy to establish) properties of the density are the 
following: 

1. An arithmetic sequence with more than one element has positive density. 

2. A sequence which grows in absolute value exponentially (in both direc- 
tions) has density 0. 

3. Z has density 1. 

4. The density of the union of two sequences (reordered to be monotonic) is 
no more than the sum of densities of the two sequences. 



Assume that 

2n 



u = 



^ diVi (3.8) 



Then in particular, for the j th coordinate vj = YluLi a i v J- We know that the 
values of <j) in the slices " higher" then those represented by u can be calculate by 
repeated iterations of the matrix A. Thus, if the last ("lower") n coordinates of 
u represent slice number 0, and the first n coordinates represent slice number 1, 
then cj) takes in slice number k the values which appear in the first n coordinates 
of Uk '.= A k u (or cquivalently - those in the last n coordinates of A k+1 u). Thus, 
for any eigenvector vl for eigenvalue A, 7^ 1, its coefficient in the representation 
of uu is Aj Oj. 



Let Xi > 1 be a maximal eigenvalue for which Oj 7^ 0. There must be at least 
one such Xi , due to the assumption of the lemma, yet - there might be several 
equal eigenvalues. Denote by t 3 the sum Xoia =a } aiVl " '> wnere w f denotes the 
j-th coordinate of vi . This may be considered as the sum of contributions of the 
eigenvectors with eigenvalue Xi to u in the j coordinate. As the eigenvectors 
are linearly independent, and a, 7^ 0, not all the t 3 's can be 0. Without loss of 
generality t 1 7^ 0. 

Note that there must be some other eigenvalue Af, < 1 such that the sum of 
contributions of eigenvectors with that eigenvalue Xi to u in the I s * coordinate 
7' 1 := X^iiA =a > a i v i l is non-zero as well. The reason is that if we consider lower 
slices (applying A ) then the contributions of eigenvectors with eigenvalues 
smaller than 1 decay exponentially and the contributions of eigenvectors and 
extended eigenvectors of 1 to the first coordinate are summed to be an arithmetic 
progression (that need not to be, at least naively, an arithmetic progression of 
integers). If there are no other elements added, then the I s * coordinate in 
"negative" slices are only the sum of such sequences. But it is very easy to 
verify that such a sum cannot be always integer (as it must be) if the exponential 
decaying sequences do not add to 0. And in our case they cannot. 

Thus, for a large enough positive integer k Uk 1 ~ t l Xi ', and for small enough 
negative integer k Uk 1 ~ r x Xb (assuming that A& is the smallest eigenvalue 
with nonnegative contribution of its eigenvectors to the I s * coordinate). Hence, 
from property 2 of densities, we have that the values that the first coordinate 
receives along its line is a set of density in Z. This analysis can be carried to 
any coordinate with nonzero contribution related to some eigenvalue other than 
1. 

Assume that there are k coordinates without any contribution from eigen- 
value A 7^ 1, and the rest n—k coordinates are contributed from such eigenvalues. 
From property 4 of the density it follows that the density of the union of these 
n — k lines is 0. On the other hand, in the other k lines, the coordinates (val- 
ues of </>) are non constant arithmetic progressions (due to the Jordan extended 
eigenvectors). Thus - they have a density (property 1). The sets of values of 
<f> in each line are pairwise disjoint (0 is bijective) and their union cover Z (cf> 
is surjective). Thus, due to properties 4,3 the density of the union of the k 
arithmetic progressions is 1. We shall now prove the following lemma, and use 
it to finish the argument. 

Lemma 3.5. If S 1 , ...S n are arithmetic sequences s.t. the density of their union 
is 1 then they cover all Z. 

Proof. Assume that S l has the form {ai + mdi} m . Let D denote the l.c.m. 
of {di}. Consider the intersection of {1,2, ...,D} with S, the union of the se- 
quences. If this intersection is all of {1, 2, ...,£)}, then S contains all of Z. 
Indeed, if k G S then k + mD is in S as well, for any integer m, as k must be- 
long to some S l , hence k + Idi must be in S for any /, but since D is a multiple 
of di, k + mD is in S as well, for any m. Thus 5 3 1, and hence S = Z. 

On the other hand, if there is some 1 < k < D not in S, a similar argument 
shows that k + mD is never in S, for any m (as the analysis above shows that 
if k + mD is in S, so is (fc + mD) — mD). But from here it is easy to verify that 
the density of S is at most 1 — -L . A contradiction. □ 



We can now finish the proof of Lemma 13.41 It follows from the lemma that 
the values of <\> in the relevant k lines (those without exponential growth) , which 
are arithmetic progressions with total density (density of the union) being 1, 
must cover Z by themselves. But then since the rest of the n — k lines also have 
values in Z, some values must repeat, in contradiction with (f> being bijective. 
And the lemma is proven. □ 

And the theorem is thus proved. □ 

Remark 1. It is easy to see that if G is just a set of points, without any edges, 
then there exists a harmonic labeling - G x Z is a collection of d copies of Z, 
name them 0, 1, 2, ..., d— 1 and then <p{a, z) := a+dz, where a £ {0, 1, 2, ..., d—1}, 
z 6 Z represent the z th point in the a copy of Z. 

Remark 2. Note that Z 2 x Z2 has a harmonic labeling and a regular tree times 
Z2 also has a harmonic labeling (easy generalization of theorems \2.2\ and \4.1\ ) ■ 
If one finds a finite graph G such that Z 2 xG doesn't have a harmonic labeling 
it also provides an example to an infinite graph G' = Z x G such that G' x Z 
doesn't have a harmonic labeling. 

Remark 3. From the proof of Theorem \S.ll one can derive a construction 
for an injective harmonic function from the ladder graph to Z. Calculating 
the Eigenvalues and Eigenvectors of the ladder's matrix A, we find there is an 
Eigenvector with an Eigenvector greater than 1 and another with an Eigenvalue 
smaller than 1 . These Eigenvectors grow exponentially in one direction of the 
ladder and decreases exponentially in the other and in opposite directions. There 
exists a linear combination of the two Eigenvectors that assures injectivity. 

4 Labeling of Z d 

Theorem 4.1. There is a harmonic labeling ofE 2 . 

Before presenting the rigorous proof we will present a sketch of the proof. 
The labeling (j> is spanned by its values on Z x {0} U Z x { — 1}. The construction 
is simply placing a rapidly increasing sequence on one copy of Z and a filling 
sequence on the other i.e. a sequence that takes care of the surjectivity of 0, 
grows slowly and jumps only over previously introduced values of <f>. It is left 
to show that <j> is injective. We get this result with tedious calculations, but the 
idea is simple. For every large element of the sequence there are only a linear 
number of other vertices with the same order of magnitude and they all have 
different coefficients in front of the smaller magnitudes. 



Remark 4. Notice that by Theorem \4- i| Z x Z admits harmonic labeling where 
as Z x [—k, k] does not admit harmonic labeling for all k (by Theorem \3.1\) . 
The essential difference between the two graphs is that a harmonic function on 
Z x [— fc, k] is spanned by a finite set and a harmonic function on Z 2 is spanned 
by an infinite set, leaving the function with more degrees of freedom that allow 
it to be bijective. See Section^ (Question^) for an open question inspired by 
this remark. 

Proof of Theorem \4-l\ We will construct a harmonic labeling of Z 2 inductively. 
To do this we need some notation. 
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Let A(x) := (x, -1) e Z 2 and B(x) := (x,0) E Z 2 . Define 

A := \jA(x) 

xez 

B := \jB(x). 



For n G IN and < i < n, set 



Set 



t/Xj, := (-n + i, i) "Upper-left" 

VR\ :=\n+l-i,i) "Upper-right" 

LL l n := (— n + i, — i — 1) "Lower- left" 

L-RJj := (n — « + 1, — i — 1) "Lower-right." 



UL(n) := |J ULi 

Q<i<n 

UR(n) := |J (7i4 

0<i<n 

LL(n) := |J LL^ 

0<i<n 

Li?(n) := |J LRi 



0<i<n 



and finally let 

S n :=( |J C/i(fc)]u[ |J L/^fc^U^ |J LL(fc)W |J ii?(fc))0 

^0<fc<n ' ^0<k<n ' ^ 0<A;<ri ' ^ 0<k<n ' 

Note that S n+1 = S„ U f/i(n + 1) U f/i?(n + 1) U LL{n + 1) U Li?(n + 1). 

Observe that, for any harmonic function / := Z 2 — > R, the value of / at the 
point (s,y) is an integer combination of the numbers f(x,y — l),f(x — l,y — 
1), /(ar+1, y—1) and f(x, y—2). Inductively, an arbitrary function / : AlJB — > Z 
extends uniquely to a harmonic function / : Z 2 — > Z. Moreover, the values of / 
on the set {A{x) : — n < x < n + 1} U {B(x) : —n < x < n + 1} determine the 
values of / on the set S n (but does not determine the value of / at any point 
inZ 2 \S„). 

We will construct a sequence of functions 

f n : {A(x) : -n < x < n + 1} U {B(x) : -n < x < n + 1} -> Z 

so that for each n, f n +i extends the function /„ and will denote by f n : S n — > Z 
the (unique) function determined by /„ using the harmonic property. We will 
do this in such a way that the common extension / : Z 2 —> Z of all of the 
functions f n is harmonic and bijective. 



The above terminology for UL(n),UR(n), LL(n) and LR(n) as "upper-left," "upper- 
right," etc. refers to its location on the boundary of the set S n . 
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Step 1 (Base of Induction): Begin by setting 

UULl) := 10! 
h{UL\) := 10!! 



h(UR° ) := 10!!! 
h(UI#) := 10!!!! 



fi(LL° ) := 

h{LL\) := -1 

fi(LR° ) := 1 

MLR ,) := 2. 

By the harmonic property, 

h(UL\) = 4-10!- 10!!- 10!!! 
fi(UR\) = 4 -10!!! -10!- 10!!!! -1 
h(LL\) = -10! 
/i(iflj) = 4 -10!!! -2. 
Step 2 (Induction Step): Suppose we have constructed the function 

/„ : {A(x) : -n < x < n + 1} U {B(x) : -n <x < n + 1} -$• Z 

so that its harmonic extension /„ : $„ — > Z satisfies (note the similarity of 
conditions 2-5), 

1. The range of f n includes each integer in the set [— n, n + 1]; 

2. For < i < n, we can estimate f n (UL l n ) by 

i 
and 

A(^ = (-l) B (/n(^-4n-/»^ 1 )+W^))±^ T /»Pti)i 

n + 1 

3. For < i < n, 

i 



f n {ULi) = {-iy{f n {ULl) - U ■ UUL^)) ± ^-f n {UL° n _S 



UiURl) = (-iy{f n (UR° n ) - 4z • AC^ij".!)) ± -— UUR ^) 

t + 1 



and 

f n (UR^^(-ir{f n (UR^-4n-f n (URl_ 1 )+f n (ULl)))±^—f n (UR^_ 1 ); 

n + 1 

4. For 1 < i < 77,, 

/„(Li;) = (-l)' t+1 (/„(^-!) - 4(» - 1) • /n(^°_ a )) ± T^-r/„(t/i°_ 2 ) 

and for n > 1, 

/n(i^) = (-l) B (/n(^^i)-4n./ n (CTL°_ 2 )+/ n (J7i^_ 1 )))±^/ n (^_ a ); 



-By a = b±c, we mean 6 — c < a < b + c 
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5. For 1 < i < n, 

f n (LRi l ) = (-lf+ 1 (f n (UR n _ 1 )-A(i-l).f n (UR n _ 2 ))± i ^ l f n (UR° n _ 2 ) 
and for n > 1, 

/ n (iK) = (-l)"(/ n (C/^_ 1 )-4n-/n([/i?°_ 2 )+/„(C/i°_ 1 )))±^ T /„(t/£°_ 2 ); 

n + 1 

6. /„ is injective. 

We will construct f n +i (extending /„) so that its harmonic extension f n +i also 
satisfies the above properties. We proceed in four stages, analyzing f n +i on the 
upper-left, upper-right, lower-left and lower-right edges of its domain, S n +\. 

Let JVi be the largest negative integer not attained by /„ and N 2 be the 
smallest positive integer not attained by /„. Let N3 be the sum of the absolute 
values of all integers in the range of f n . We define 

f n+1 (LL° n+1 ) := m 

fn+i(LR n+1 ) := N 2 

f n+ i(UL° n+1 ) := lO»>-»r+»* 

f n+1 (UR° n+1 ) := 10 10 - 

and fn+i = fn for all points where /„ is defined. Therefore f n +i = fn on S n . 
By construction (of / n +i), property one is satisfied for f n +i- 

Step 2.1 (Upper-left edge): We show that property 2 is satisfied by induction 
(on i). Note that f n +i satisfies property two at UL l n+l because 

f n+1 (Ul} n+1 ) = 4/„ +1 (C/L, 1 l ) - f n+1 (UL° n+1 ) - U+^UL ^) - / n+1 (iL°_ 1 ). 

Now suppose 1 < i < n and f n +i satisfies property two at UL 3 n+1 for each 
< j < %■ Using the harmonic property, 

f n+1 (UL^\\) = 4f n+ i(UL l n ) - f n+l {UU n+l ) - f n+1 (UL l n ^) - f n+1 (UV-_\). 

(4.1) 
By property two and the induction hypothesis, 

f n+l {UU n ) = (-iy{f n+1 (ULl) - U ■ f n+1 (UL° n _ t )) ± ^J-f^ULl,!) 

' % + 1 

f n+1 (ULi l+1 ) = (-iy(f n+1 (UL° n+1 ) - 4* • f n+ i{UL a n )) ± -l-f n+1 (UL Q n ). 

f (TIT ® } 

By construction of /„+i(C/L°_ 1 ), Z +1 frjfjrx < 10 _ ". Therefore (recalling that 
i < n), 

-/n+i(17i£_i) - f n+1 (ULl-_\) - (-l)Mi • f n+1 {UL a n ^) ± l —}-f n+1 (ULl^) 
<(4n + 5)f n+1 (UL° n _ 1 ) 

(4.2) 
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< -^—fn+l(ULl) 

n z + 1 



Combining this with (|4.ip shows that property two holds at UL l ^ x . By induc- 
tion, property two holds for < i < n. A similar argument shows the case i = n. 

Step 2.2 (Upper-right edge): Property three follows by a similar argument to 
property two using UR l n+l instead of UL\ 1+l . 

Step 2.3 (Lower-left edge): Note that if n > 1, 

fn+i(LL n ) = 4f n+ i(LL n _ 1 ) - f n+ i(LL n ) - f n+ i(LL n _ 2 ) - f n+1 (UL n _ 1 ) 

(4.3) 

= -f n - 1 (UL° l _ 1 )±3n 2 . (4.4) 

The rest of property four follows by a similar argument to property two. 

Step 2.4 (Lower-right edge): Property five follows from a similar argument to 
property four. 

Step 2.5 (Injectivity) : Recall that f n +i is injective when restricted to S n . We 
need only see that it remains injective on its whole domain S n +i. Property two 
shows that f n+ \ is injective when restricted to UL(n + 1) (and properties 3-5 
show that it is injective when restricted to UR(n + l), LL(n + l) and LR(n + l), 
respectively). 

If (x, y) G UL(n + 1), then property two shows that 

\f n+1 (x,y)\ > \f n+1 (UL° n+1 ) - An ■ f n+1 {UL° n )\ > 2- max |/ n+ i(o,6)| 

by construction of f n+1 {UL° n+1 ) (which is io^(a, b) esJf"+i( a > b )\y Therefore f n+1 
is injective when restricted to S n U UL(n + 1). Similarly it is injective when 
restricted to S n U UR(n + 1), S n U LL(n + 1) and S n U LR(n + 1). 
We claim that f n +i is injective on all of SVi+i- By properties 2-5, 

If (x,y) e UL(n + 1) then f n+1 (UL° n+1 ) = 2 ±J f n +i(UL° n+1 )± 

If (x,y) G UR(n + 1) then f n+1 (UR° n+1 ) = 2 ±1 f n+1 (URl +1 ); 

Ii(x,y)ELL(n + l) then f n+1 (LL° n+1 ) = 2 ±1 f n+1 (UL° n ); 

If (x,y) E LR(n + 1) then f n+1 (LR° n+1 ) = 2 ±1 f n+1 (UR° n ). 

Recall that, by construction, 

f n+1 (LRl +1 ) > lo/^C^iUi) > io^" +1<c/R " ) > lol° lo/ " +l(C/i ^ 

Thus, fn+i is injective on S n +i (i.e., property six holds). 

Step 3 (Construction of labeling): By induction, we can construct the functions 
/„ and /„ for all n so that they have the six properties above. Since f n +\ 
extends /„ we can dehne f(x,y) to be f n (x,y) for any n large enough that 
(x,y) G S n - By construction, f n +\ is harmonic on S n so, by properties one and 
six, / is a harmonic labeling. □ 

Theorem 4.2. There is a harmonic labeling for ~L d for all d G N. 



11 



Proof sketch. After reading the proof of theorem |4. II one can see how to general- 
ize the proof for d > 2 but at the same time it seems that writing a rigorous proof 
is rather complicated and tedious. The labeling of Z d is spanned by its values 
on Z d_1 x {0} U Z rf_1 x { — 1}. We construct 4> on the boundary of [— n,n] d ~ 1 . 
For d = 2 we added the rapidly increasing sequence in two points each step 
(the boundary of [— n, n]), for d — 3 wc introduce the sequence to the boundary 
of a rectangle and so on. Lets take d — 3 as an example of how to continue. 
By introducing the rapid sequence on a rectangle boundary of one slice and a 
filling sequence on another slice we can extend to a pyramid shape by using 
the harmonic property. Now Each order of magnitude dominates a face of the 
pyramid, where again the equations are different for each vertex. Proving the 
difference of equations is the technical part to rigorize, but it is not hard for one 
to convince himself the construction works for all dimensions. The trick is to 
take the increasing sequence blow fast enough such that the increase will be fast 
enough to beat the polynomial number of elements of the lower magnitudes. □ 

5 Open problems 

1. Geometric properties of label spanning set. We say that a set V C V(G) 
is label spanning if using the harmonic property, given any values of <p 
on V we can extend <p to the entire graph. We saw several examples of 
finite label spanning sets (Cross, Z and the ladder). The first interesting 
question is are there any graphs other then several disconnected copies of 
Z which admit harmonic labeling spanned by a finite set? What can one 
tell about the graph from its spanning set? 

2. Growth rate of harmonic labeling. Set some vertex as the origin 0. Look 
at the set A n = {<f){x) : d(x, 0) < n}, where d is the graph distance. What 
can one learn about a graph that admits harmonic labeling by the growth 
rate of A n l 

3. Are there any infinite graphs G such that G x Z doesn't admit harmonic 
labeling. Does any G x Z where G is infinite transitive admit harmonic 
labeling? 

4. Is there a Cayley graph, which is not a finite extension of Z, that doesn't 
admit harmonic labeling? What about the Hcisenberg group? Other 
natural graphs to consider are hyperbolic lattices. 

5. In Remark [3] we saw an example of a graph that does not admit harmonic 
labeling but has an injective harmonic function to Z. Does any vertex 
transitive graph admit an injective harmonic function to Z? Note that 
Dropping the injectivity constraint is easier as every graph that doesn't 
have a finite label spanning set admits a surjectivc harmonic function to 
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